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ABSTRACT 

Multi-wavelength observations are mandatory to understand the physical 
properties of astrophysical sources. In this paper we use observations in the far 
infrared to radio range to derive the spectral energy distribution (SED) of the 
Carina nebula. To do this, we carefully subtract the irregularly varying diffuse 
emission from the Galactic plane, which can be of the order of 10 % of the nebula 
flux at these wavelengths. We find that the far infrared SED can be modeled as 
emission from a dust population with a single temperature T d = (34.5^^g), and 
with a spectral index of emissivity a = — 1.37loo8- We also find a total infrared 
luminosity of the Nebula of (7.4l^ 4 ) x 1O 6 L and, assuming a single temperature 
of the dust, a mass of the dust of (95OOl|^)M . 

Subject headings: Infrared: ISM — Instrumentation: photometers — Methods: data 
analysis — Methods: statistical — radio continuum: ISM — Submillimeter: ISM 
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Introduction 



The Carina nebula (NGC 3372 or RCW 53) has one of the highest surface brightnesses 
of any nebula in the southern hemisphere. It contains the most luminous and one of the 
most massive stars of our Galaxy: eta Carinae, a LBV star. Its distanc e (2.3 kpc) is known 



accurat ely, mainly t 



nebula ((Smith 



19931 ) of the Homunculus 



; irough the expansion parallax (lAUen fc Hillierl . 
20051 ). The nebula has a particularly low extinction with respect to other 
star formation regions, and presents an anomalous reddening law ( Smith Il987 ). The 
nebula hosts, moreover, more than 60 O stars, three WNL stars, an 02 s upergiant, an d 



a B1.5 supergiant which suggests the presence of a proto OB association ( (Smith 



20061 ). 



The Carina nebula offers a unique laboratory to study different astrophysical problems: 
star formation, feedback on the star formation provided by young massive stars, as well as 
triggering of a second generation of star formation from the most massive ones; small-scale 
phenomena (evap orating protoplane tary disks, jets, cometary clouds, dust pillars, globules, 



Herbig-Haro jets 



Sm ith et al. 



molecular clouds (IRathborne et al. 



200 41), an d photodissociation regions on the surfaces of 



20021 )). High-quality data spanning from infrared (IR) 



to radio wavelengths show a heavily active star formation region in which the LBV star, as 
well as other massive ones, play a crucial role. 

Knowledge of the nebula improves our knowledge of the interstellar medium, 
and, in particular, our models of interstellar d ust. Observational cosmology has an 



interest in understanding dust properties too (IDeschenes 



2008 



Lazarian fc Finkbeiner 



2003 



Masi et al. 



Dunkley et al. 



2008 



20011 ): current measurements of the polarization 



of the cosmic microwave background (CMB) are limited by our limited ability to 
subtract foreground emission. Future experi ments like the Planck and Hersc hel satellites 



(launched on 2009 May 14), and t 



2004 



Reichborn-Kjennerud et al. 



ie PILOT ([Bernard et al 



2010), SPIDER flMontrov et al. 



2007a 



b|), EBEX f Oxlev et al. 



20061 ). and Balloon 
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Observations of Millimetric Extragalactic Radiation and Geophysics (BOOMERanG) 
balloons, with their polarization surveys of the sky will provide data to understand dust 



properties: like the emissivities and t 



polarized emission ( iHildebrand et ai- 



re wavelength-dependence of the far infrared (FIR) 



1999 



Lazarian fc Hoang 



2007: 



Martin 



20071 ). This 



knowledge is very much needed for a better extraction of the finest CMB polarization 
properties, which can provide information on cosmol ogical inflation, in post- Planck 



missions like CMB-Pol ( iBaumann et al. 



20091 ). B-Pol fide Bernardis et al. 



2009|) and COrE 



(http: / /www. core-mission. org). 

The paper is organized as follows: in Section^ we briefly summarize the current 
knowledge of the Carina nebula. In Section[3] we present the data we used for our analysis. 
The subtraction of the background is shown to be, with the calibration error, the dominant 
source of error for the measured fluxes (SectionHJ). In Section^ we present the FIR to radio 
spectral energy distribution (SED), and the best fit we performed using y 2 minimization 
and a Monte Carlo Markov Chain Metropolis-Hastings (MCMCMH) algorithm for the 
estimate of the confidence intervals. We derive some global properties of the nebula and 
compare our results to previous literature data in Section[6j and we conclude summarizing 
the results and future activities in Section[0 



2. The SED and the global properties 



Smith fc Brooks! (120071 ) have integrated IRAS (FIR) data and the Parkes (radio) data 
to estimate the flux of the Carina nebula in a box 1°58' x 2°45' wide. The resulting SED 
is described by three graybody components having temperatures of 220, 80, and 35 K with 
emissivity oc A -1 , and by an optically thin thermal bremsstrahlung radio continuum emission 
component (oc v~ 0A ). Their analysis used the two data at 60 and 100 /im to fit the Carina 
nebula SED in the wavelength range (20-5000) /mi, with two free parameters (temperature 
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and emissivity) , assuming a spectral index of emissivity -1; the single measurement at radio 
wavelengths is used to measure the amplitude of the thermal bremsstrahlung component. 

Luminosity at these wavelengths is independent of dust grains radiu s and emissivity 



1984T ). With a typical 



as long as the radius of the grains is small, < 0.2/im ( iDraine &: Led 
grain density, p, of about 3gcm~ 3 and a temperature T, the dust mass which emits an IR 



luminosity, Lm, is ( jSmith fc Gehrz 



20051 ): 



M dust {M & ) 



122p 



3a SB TS 

with a sb being the Stefan-Boltzmann con stant expressed in cgs units. For the dust 



(1) 



population at 35 K 



Smith fc Brooks! (120071 ) found a luminosity and a mass equal to 



7.6xl0 6 L^ and 9600 



by 



Mq. respectively. In the following we keep the same area considered 



Smith fc Brooksl (120071 ) and we estimate the integrated flux of additional experiments 



carried out at different wavelengths. 



Studies of the Carina nebula at wavelength s shorter than 60/xm ar e 



1995 



Morris et al. 



1999 



Rathborne et al. 



2002 



already present 



Smith fc Brooks 



in lit erature (I Cox et al. 
20071 ). In this work we demonstrate that in the FIR portion of the spectrum the B98 



and Wilkinson Microwave Anisotropy Probe ( WMAP) data improve our knowledge of the 
nebula. 



3. Multiwavelength observations 

In this section, we describe the multiwavelength data we used and how they have been 
analyzed. Although we used different experiments, some parts in the data analysis are the 
same. 

For each experiment the average value of the background has been estimated as 
described in SectionHl For the background estimation, we have also considered a possible 
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contribution from CMB anisotropy, but we found that it is negligible, given the large 
number of independent sky regions present in our integration box. 



3.1. COBE 



The DIRBE instrument of the COBE satellite has produced infrared absolute sky 
brightness maps of the Carina nebula i n the wavelength ra nge 1.25 to 240 /im with an 



instantaneous field of view of0.°7x0.°7 ( IHauser et al 



1998al ). For our analysis we used data 



at wavelengths longer than 60/im. The quadrilateralized spherical cube projection of the 
DIRBE ma ps has been convert ed into the Hierarchical Equal Area isoLatitude PIXelization, 



HEALPIX ( lOorski et al. 



20051); this has required a change of pixel size with respect to the 



original: from about O.°32o to 27/5 (nside=128). From these d ata we have calcu l ated the 



maps and the integrated fluxes in the same box considered by iSmith fc Brooks! (120071 ) . 
The results are shown in Figured] and TableHl The errors in the calculated fluxes are 
the quadrature sum of the uncertainty in the absolute gain calibration, the uncertainty 
in the instrumental offset (zero point), the standard deviation of the weighted average 



photometry values ( IHauser et al 



1998bl ). and the contribution from both the background 



diffuse emission of the Galaxy (see Section!!]) and the rms CMB anisotropy. 



3.2. IRAS 



IRAS has observed, at wavelengths of 60 and 100 /im, the Carina nebula with an 



angular resolution of 1/5 x 4/7, and 3/0x5/0, res pectively (INeugebauer et al. 



have used the latest processed IRAS images IRIS (IDeschenes &: Lagache 



19M)- We 



20051 ): they have 



an absolute calibration coherent with DIRBE, a bet ter zodiacal light subtr action, and better 



destriping with respect to the original ISSA maps ( iWheelock et al. 



19931 ) . A preliminary 
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comparison with the original ISSA data gave flux values 22% and 7.9% smaller than the 
IRIS data at 60 and 100 /im, respectively, the discrepancy with the IRIS data is 1.5 (0.3) a 
at 60 (100) /im. For coherence with the other experiments, that are i n the HEALP IX 



pixelization, we used the IRIS data. The resolution is 1/7 (nside 2048) (IGorski et al.ll2005l ). 
The brightness map in our integration box is reported in Figured] and the resulting flux in 
TablegJ 



3.3. BOOMERanG 98 



The stratospheric balloon experiment BOOMERanG 98, hereafter B98, (ICrill et al. 



20031) observed the Carina nebula at frequencies of 90, 150, 220, and 410 GHz ( iCoble et al. 



20051 ) with an angular resolution of about 18/1, 1 1/7, 15/9 and 15 /0, respectively. For 



the 410 GHz chann el the calibration error is 25% (IMasi et al. 



flRuhl et al. 



200 lh . for the others 10% 



20031 ). Since the B98 survey does not fill completely the Smith integration 



box, the missing region (white region in B98 maps, Figured]) is filled with the average signal 
from an observed background region of area 0.°2x0.°2 (see Tab. [I]). The missing area is at 
most 10% of the integration box, in a low background region. From the WMAP7 map at 
a similar frequency (94 GHz) we see that neglecting completely the contribution from this 
area the flux decreases by 1.2%, so the error we introduce with our procedure is certainly 
less than 1% and in any case is smaller than the errors introduced by the background 
subtraction procedure (see SectionH]). Fluxes have been corrected for the instrumental 
response of the telescope using the transmissions functions of the 16 B98 filters. The signal, 
in Jansky, for each channel i, Si, is given by 

kA,iTcMB 



S; 



Le^B^^AQdu J A 



ei{y)Ii{y)dv, 



(2) 



in which Tqmb = 2.73 K is the temperature of the CMB (IFixsen et al. 



19961 ). x = hv/(kT) 



(with k being the Boltzmann constant and h being the Planck constant), kA,i is the 
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calibration factor of the ith channel, eiiv) the filter transmission function of the zth channel, 
and l{{y) is the theoretical emission model. 

Table 1: Regions Observed by B98 in Different Frequency Channels, and Used to Fill Miss- 
ing Regions in the B98 Survey; the Columns List Respectively the B98 Channels, and the 



Coordinates of the Edges of the Regions: t 


mini ""max-. 








Channel 


*"min 


I 


max 


bmin 


bmax 


B90A, B90B 


287°.9 


21 


S8M 


-1°. 


.8 


-1°.6 


B150A, B150B, B150A2, B150B1, B150B2 


288°.3 


21 


38°. 5 


-2°. 


1 


-r.9 


B150A1 


287°.6 


21 


37°. 8 


-2°. 


.1 


-1°.9 


B220A1, B220B1, B220B2 


287°.9 


21 


38°. 1 


-1°. 


.8 


-1°.6 


B220A2 


288°.3 


21 


38°. 5 


-2°. 


.1 


-1°.9 


B410A2, B410B1 


288°.3 


21 


38°. 5 


-2°. 


1 


-1°.9 


B410A1 


288°.3 


21 


38°. 5 


-1°. 


.8 


-1°.6 



We have assumed a typical power law, Ii{v) = Io,i{v / Vo,iY , where Io,ith is the flux 
measured at frequency uq^, the central frequency of the ith channel. The integration is 
performed over the entire band of the filter, Avi. The color corrections (TableEJ) are less 
than 1% for an optically thin thermal bremsstrahlung (/3 = —0.1), Rayleigh- Jeans spectrum 
(/3 = 2) or a power law with f3 = 3. Since in the B98 wavelength range (excluding the 
90 GHz channels) we expect thermal dust emission to be predominant over bremsstrahlung 
emission, we carried our analysis with the color corrections relative to = 2. B98 was a 
scanning instrument. To remove instrumental offsets and drifts, we have high-pass filtered 
the time-ordered data. The cutoff of the high-pass filter corresponds to an angle 9c in the 
sky, so the B98 maps we have used do not contain structures at scales larger than 6 C - We 
have verified that the fluxes in our target region are not significantly affected by filtering if 
6 C > 10° (see Tabled . 
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Table 2: Variations of the Integrated Fluxes (Color-corrected) Measured at 150 GHz for 
Different Spectral Index of SED, /3, and Different Cutoff $c of the High-pass Filter Applied 
to the Time-ordered Data. The fluxes are corrected for the background emission (n = 5, see 
Sectioni]). 



Filter 





Flux (Jy) 


5°.0 


-0.1 


920 ± 150 




2.0 


920 ± 150 




3.0 


910 ± 150 


10°.0 


-0.1 


1630 ± 190 




2.0 


1620 ± 190 




3.0 


1610 ± 190 


15°.0 


-0.1 


1920 ±210 




2.0 


1920 ±210 




3.0 


1900 ±210 
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The errors in the fluxes are the quadrature sum of the standard deviations of the 
measured signals, the calibration uncertainties, and the Galactic background contributions. 
Depending on the channel, the calibration error dominates on the background one and vice 
versa (TableHJ) while the statistical error represents a small contribution to the total. 

The final integrated fluxes are computed as averages of different channels with the same 
frequency: one for 90 GHz, four for 150 GHz, three for 220 GHz, and three for 410 GHz. The 
results are reported in TableH] and Figured! 



3.4. WMAP 



20031 ) at 



The Carina nebula has been observed by the WMAP satellite (IBennett et al. 
the frequencies: 23 GHz (if-band), 33 GHz (Ka-band), 41 GHz (Q-band), 61 GHz (F-band) 
and 94 GHz (VF-band); the seven-year release data, hereafter WMAP7, provide maps of 
the nebula w ith FWHM angular r esolu tion of about 0°.81, 0°.62, 0°.48, 0.°33, and 0.°21, 



respectively (IHinshaw et al. 



2007 . 



20091 ) . We used the temperature (Stokes I parameter) 



maps per frequency band obtained from the Lambda archive ( |http: / /lamb da. gsfc.nasa.gov/ ). 
The fluxes have been corrected for the instrumental responses of each of the 20 differential 
radiometers; we averaged the bandpass frequency responses of all the radiometers in each 
frequency band. The errors have been estimated as the quadrat ure sum of statistical error, 



background uncertainty, calibration error (0.5%) and pixel noise (IHinshaw et al. 



2009|). The 



results for the integrated fluxes are reported in TableH] and the maps in Figure[TJ 



3.5. Spitzer 



The MIPS instrument on board of the Spitzer satellite (Spitzer archive: 
http://archive.spitzer.caltech.edu/) has observed the Carina nebula at far infrared 



wavelengths of 24, 70, and 160 /im (IRieke et al. 1 120041 ) - However, the photoconductors of 
MIPS were significantly saturated by the strong emission of the nebula, so we decided to 
exclude these data from our analysis. 



3.6. Multiwavelength maps 

The maps resulting from our analysis of all the experiments are displayed in Figured] 
they are labeled with their HEALPIX pixel scale. The COBE maps, nside=128, have been 
filtered with a Gaussian low-pass filter (FWHM=0°.35). Because of the low resolution of the 
instruments, the inner morphology of the Carina nebula is not resolved in the COBE-DIRBE 

'he lower resolution of the 



Smith fc Brooksl (120071 ). does 



maps and in the three lower frequency channels of WMAP. 
B98 instrument, with respect to the experiments analyzed by 
not allow a real pixel to pixel comparison: there is, nevertheless, significant morphological 
resemblance between the IRAS-IRIS, the B98 maps and the 61 and 94 GHz channels of 
WMAP7. 



4. Background subtraction 



The region around the Carina nebula, near the Sagittarius-Carina spiral arm, is rich 
in massive star formation and HII regions: NGC 3603, NGC 3576, NGC 3572, NGC 3324, 



NGC 3293 and RCW 49 flSmith fc Brooks 



20081 ). With respect to these regions Carina is 



larger, brighter, closer (2.3 kpc against 6-8 kpc) and has a lower extinction and reddening; 
therefore, its emission should dominate. 

Background contributions (CMB, sky and irregularly varying diffuse emission) have 
been studied averaging flux values with two different numbers of background regions 
(labeled with the index n). In both cases we use six sets of reference regions, each set 
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including five regions (n = 5) and two regions (n = 2). All regions have an area of 0.°5 x 0.°5 
(see Tabled]). For a given n, the reference regions are labeled with the index m. For n = 5 
the regions have been randomly chosen around the Carina nebula, for n = 2 the regions 
are close to the nebula and symmetrically located with respect to its center. Only regions 
observed by all the experiments and close to the Carina nebula have been selected. We 
computed the average of the fluxes of all the regions in each set. We used the average of all 
the sets as an estimate of the background, and its standard deviation as an estimate of the 
error on the background. 

Table 3: Selected Sky Areas for the Study of the Diffuse Background Contribution. Top 
part: n=5, bottom part: n=2; m is the index of the region, bi (2=1,5) are the galactic 
coordinates (in degrees) for the centers of the regions. 



m 




ti 


h 




£2 


b 2 




4 


h 




h 




h 


h 


1 


21 


39.75 


1.25 


21 


35.15 


-3.75 


21 


35.15 


1.25 


287.20 


1.25 


21 


35.75 


-2.35 


2 


21 


39.75 


2.25 


2 C ' 


36.15 


-3.75 


2< l 


35.15 


2.25 


287.20 


2.25 


21 


35.75 


-3.35 


3 


21 


39.75 


3.25 


21 


34.75 


-2.25 


2< l 


38.25 


1.25 


284.25 


-3.75 


21 


35.55 


3.25 


4 


21 


38.75 


1.75 


21 


36.25 


1.75 


2< l 


34.25 


-4.25 


285.25 


-2.25 


21 


37.75 


3.75 


5 


21 


38.25 


3.25 


21 


34.25 


-2.75 


21 


S4.25 


1.75 


289.75 


4.75 


21 


39.25 


1.75 


6 


21 


36.25 


4.75 


21 


34.75 


-3.25 


2! 


34.25 


1.25 


290.25 


2.75 


21 


34.75 


2.95 












m 


ii 






h 


b 2 




















1 


285.25 




-3.75 


289.75 


2.25 




















2 


286.25 




-2.25 


288.75 


0.75 




















3 


284.75 




-3.25 


290.75 


1.75 




















4 


284.25 




-3.75 


290.75 


2.25 




















5 


285.75 




-3.25 


289.25 


1.75 




















6 


285.25 




-2.75 


289.75 


1.25 
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An indicator of the effect of background fluctuations is the ratio between the standard 
deviation of the average background and the nebula flux (corrected for the background). 
When computed for all the regions of Tabled this indicator has a trend increasing with 
wavelength, and is maximum between the 240 jttm DIRBE and the 410 GHz B98 images; on 
the contrary this indicator is basically independent of the considered sky region (FigureE]). 

The background contribution is more dependent on the number of regions considered 
for the analysis; this seems to be true for all the experiments. Similar background trends 
are found considering a small number of sky regions (maximum three) with different solid 
angles (from 1° x 1° up to 6° x 2°). 

In both cases, the WMAP data are the least affected by the background, followed 
by B98, IRAS and COBE. The background disagreement between the B98 90 GHz and 
WMAP5 94 GHz is due to the high-pass filtering of the BOOMERanG data. 

Our analysis seems to indicate, for the irregularly varying emission surrounding the 
Carina nebula, a significative spectral dependence and a large scale nature (larger than 
6° x 4°). 

The error on the background in case n = 5 is about 1.5 times smaller than in case n=2. 
Since case n = 5 has better statistical significance, we have selected this configuration for 
the subtraction of the background in our integrated flux analysis. However, we have verified 
that all the results we give below remain very similar if we use the background estimates of 
case n = 2. 

5. SED of the Carina Nebula, from FIR to cm- wavelengths 

The integrated fluxes, with the background subtraction described in SectionHl are 
reported in TableSJ In Figure|4]we plot the SED of the Carina nebula. At 60 and 100 /zm, 
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Table 4: Integrated Fluxes for n=5. For a given wavelength of observation, A, of an experi- 
ment, we report, in the order: experiment, flux, total, statistical, background, and calibration 
errors. 



A 


Experiment 


Flux (Jy) 


Error (%) 


Stat (%) 


Bckg (%) 


Cal (%) 


60 jum 


IRAS (IRIS) 


(7.5 ±0.8) x 10 5 


10.4 


0.04 


1.0 


10 


60 /xm 


COBE 


(5.7 ±0.6) x 10 5 


10.5 


0.2 


1.4 


10 


100 fim 


IRAS (IRIS) 


(9.7 ± 1.3) x 10 5 


13.7 


0.03 


2.2 


13 


100 fjtm 


COBE 


(8.0 ± 1.1) x 10 5 


13.8 


0.2 


2.6 


13 


140 /iin 


COBE 


(9.0 ± 1.0) x 10 5 


11.4 


0.1 


4.1 


11 


240 jum 


COBE 


(4.1 ±0.5) x 10 5 


13.1 


0.1 


6.1 


12 


731 /iin 


B98 


(1.6 ±0.4) x 10 4 


26.6 


0.2 


9.2 


25 


1.4 mm 


B98 


(4.2 ±0.5) x 10 3 


12.8 


0.3 


8.0 


10 


2.0 mm 


B98 


(1.6 ±0.2) x 10 3 


11.8 


0.6 


6.3 


10 


3.2 mm 


WMAP7 


(10.2 ±0.2) x 10 2 


2.4 


1.1-10 3 


2.3 


1.7 


3.3 mm 


B98 


(0.9 ±0.1) x 10 3 


11.8 


0.5 


6.2 


10 


5.0 mm 


WMAP7 


(9.5 ±0.2) x 10 2 


1.9 


9.7-10- 4 


1.8 


1.6 


7.3 mm 


WMAP7 


(9.6 ±0.2) x 10 2 


1.9 


5.M0~ 4 


1.8 


1.6 


9.7 mm 


WMAP7 


(9.9 ±0.2) x 10 2 


1.9 


5.4-10" 4 


1.9 


1.6 


1.3 cm 


WMAP7 


(9.7 ±0.2) x 10 2 


2.0 


2.6-10" 4 


2.1 


1.6 
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the fluxes from COBE and from the Z/L45'-IRIS are mutuall y consistent with i n 2 an d 1 a, 
respectively. The JiMS'-IRIS flux at 60 /im is la away from ISmith &: Brooks! (120071 ). and 



the COBE one is at 2.9a. The IR AS -IRIS and the COBE fluxes at 100 /im are at 2.4 and 



3.5er, respectively, with respect to 



Smith fc Brooks! (120071 ). The extrapolation at 3.4 cm 



from B98 90 GHz and all WMAP7 data is ab out 10 3 Jy, at 1.6a from the Parkes data at 



the same wavelength (ISmith fc Brooks! 120071 ). This could be due to either an absolute 



calibration / offset problem in the Parkes data, or the presence of a further background 
rapidly decreasing with frequency. 

We fitted the SED, excluding the COBE points at 60 and 100 /mi, using the sum of a 
graybody (with an emissivity power law (A/Ai) ai and amplitude ei) and an optically thin 
bremsstrahlung component: 

2hu 3 1 



(3) 



where F v is the flux density, in Jy, and f2& the area, in sr, of the integration box. A 
chi-square minimization with four free parameters {x 2 'A) (the emissivity amplitude at 
Ai = 100 /im, e%; the emissivity spectral index a x ; the temperature of the graybody 
component T\, and the amplitude of the free- free component a s at the Parkes frequency 
u ) gives: X 2 reduced{x 2, 4) =1-16, e 1(x2;4) = (1.09 ± 0.13) x lO" 3 , T 1(x2;4) = (35.0 ± 2.0) K, 
ai( x 2. 4 ) = —1.32 ± 0.09, and a s ( x 2 ;4 ) = (1100 ± 21) Jy. The bremsstrahlung amplitude is well 
constrained; the IRAS / COBE data bind the peak, while the B98 ones constrain the long 
wavelength trend of the graybody law. 

We studied the correlation between parameters using an MCMCMH algorithm 



( lAndrieu et al. 



2003 



Lewis fc Bridle 



2002 



MacKay et al. 



20021 ). The resulting contour 



plots, showing Confidence Levels (CL) regions of 68.3%, 95.4, and 99.7% are reported in 
Figure|3j 



Our parameters are well determined even without the introduction of priors in the 
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algorithm. The graybody parameters are, as expected, not correlated to the bremsstrahlung 
one; the temperature, the emissivity amplitude and its spectral index, instead, are partially 
correlated. The estimates for the parameters (mcmc; 4), calculated as the average value 
of the random numbers accepted by the algorithm, are: ei( mcrnc;4 ) = (l.lOlg^g) x 10~ 3 , 
Ti(mcmc;A) = (34.5ti; g ) K, ai( mcm , c;4 ) = (— 1.37io!o8); an d a si(mcmc;4) = (llisijj) Jy with a 
reduced Xreduced(mcmc-i) =1-009. As expected the MCMCMH provides a better chi-square 
since the algorithm is not limited by the resolution of the parameters grid step as in a 
simple grid-based chi-square minimization. 

We also investigated the possibility that the SED of the Carina nebula could be 
described by two graybodies: we added the free parameters e 2 , T 2 , and a 2 , normalizing 
the second emissivity power law at A 2 = 2000 /im. The minimization of the x 2 - sc i u are fit 
provides hints for the possible existence of a colder dust component, with temperature 
of about 18K and emissivity of about 2xl0 -5 , while the warmer component would have 
a temperature of about 40 K, with a lower emissivity with respect to the case of a single 
population. 

However, the results are unstable because of the large number of free parameters 



with respect to the number of indep endent data and t 



temperature and the spectral index (jDupac et al. 



re intrinsic degeneracy between the 



20031 ); moreover, the improvement of the 



fit with respect to the four parameters case is not significant. 

We used these results as input for an MCMCMH run with the same seven free 
parameters. We restricted the parameters space, adding uniform priors for the conditional 
probabilities; therefore we set T 2 > 0, e 2 > and a 2 varying in the range [-2.4 ; -0.8]. The 
latter range has been selected following the temperature depe ndence of the submillimeter 



spectral index derived by PRONAOS data (jDupac et al. 



2003|) 



The graybody parameters still remain not correlated to a s , but the small number of 
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experimental data, with respect to the free parameters, is not sufficient to constrain the 
parameters of the second graybody and to study their correlation properties. As a check 
of the instability of this result, we performed 100 MCMCMH runs, with the same nstep 
and burn-in, and compare the results of the first run to the other 99 ones. Large variations 
were found, for all but the a s parameter, and, in any case, the lowest x 2 value we found was 
larger than the previous one. 

We tried to elim i nate the parameter a 2 using its de pendence from T 2 (see, e.g., 



Paradis et al. I ( 120101 ) 



Veneziani et al. 



(|2Qld)); following iDupac et all (120031 ) we tried 



"2 = — ■ (4) 

0.4 + 0.008T 2 v ; 

An MCMCMH with six parameters confirms the results of the MCMCMH with seven 
parameters of a warmer dust component with a lower emissivity and a second colder dust 
population with temperature of about 18 K and emissivity of about 2xl0 -5 . This is not a 
significant improvement over the four parameters fit, given the reduction of the degrees of 
freedom. 

We conclude that the data, while not excluding it, do not require an additional cold 
dust component. 

We also checked how the background removal method affects the best fit parameters. 
In case n = 2 we find a reduced \ 2 lower by 10% with respect to case n = 5, and a 
temperature 7\ warmer by 2K. All the best fit parameters are consistent with ones obtained 
in case n = 5 within la. 

Integrating the spectral flux density describing the 34.5 K dust population and 



considering an average distance of the nebula equal to 2.3 kpc (lAllen fc Hillierl 119931 ). we 
calculate the total IR luminosity for this dust population and the corresponding dust mass 
(SectionEJ) obtaining (7.4+H) x 10 6 L Q and (95OO±fjj[£)M ; both are well consistent with 
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the corresponding estimates of 



Smith fe Brooksl (120071 ). In presen ce of two dust grains 



Smith fc Brooksl (120071 ) 



populations the total luminosity is consistent, within la, with the 
estimate. On the contrary, the mass determination is strongly dependent on the assumption 
of a single dust population. Two dust populations hypothesis does not affect the luminosity 
because the SED remains very similar, but has a significant impact on the mass estimation 
as the mass strongly depends on the temperature. The nebula mass estimate increases up 
to 9.2 xlO 5 solar masses if we add a second population with temperature of about 18 K (see 
above). 



6. Discussion and future perspectives 



The determination of the SED improves our knowledge of this star forming region. In 
fact, the spectral index is a signature of a given type of dust (amorphous carbon and/or 
porous aggregates of silicates and graphite rather than silicates in ice mantles), and the 
parameters describing the SED are used to estimate the luminosity and the mass of a given 
dust population. 



Compared to ISmith fc Brooksl (120071 ). our results agree with the presence of a dust 
population with temperature of 35 K, although the fits with a double dust population might 
suggest the presence of slightly warmer dust (about 40 K instead of 35 K). This results in 
a significantly lower es timate of the mass for the warm dust component. At variance with 



Smith fc Brooksl (120071 ). in our analysis the emissivity spectral index is not as sumed a prio ri, 



but i s obtained from the data. A spectral index around -1.4 was predicted (IKoike et al. 



19801 ) from laboratory measurements on amorphous carbon, although this can be produced 



also by silicates, amorph ous carbon and graphite bounded together in a porous aggregate 



(IMathis fc Whiffen 



19891 ). A spectral index around -1 was also measured at 1.3 mm for 



observations towards two pre-main-sequence stars in the Taurus dark clouds; there thermal 
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emission is due to sm all particles embedded in disks of gas and dust surrounding the stars 



( IBeckwith et al 



19901 ) . Values between (-1 - 0) and (-2 - -1) are found, respectively, in 



submi 



limeter and millimeter wavelength observations toward T Tauri and FU Orionis 



stars (IWeintraub et al. 



1989|) ; the former range is related to extended dust clouds, the 
latter range is typical of dust in molecular clouds and compact HII regions. PRONAOS 
observations of Ori on, Ml 7, p Ophiuchi and Cygnus B show a spectral index around -1 



(IDupac et al 



2003 ), again consistent with ou r results. Early Planck results for diffuse dust 



in our Galaxy (IPlanck Collaboration 



2011al ) show a spectral index (-1.8), steeper than our 



value. However, the dust population sampled here is in a different (warmer) environment. 



In fact, indications of a steepening of the spectral index for c older temperatu r es, li 



KC 



in Equation (|H), have 



Veneziani et al. 



Deen noted by several analysis (see e.g. 



fcoioh ). 



Paradis et al. 



bold); 



For our 4 parameters fits, our best fit spectral indexes are: auy2 



(x 2 ;4) 



-1.32 and 



o>i(mcmc;i) — —1-37. Using the temperatures estimated from the fit, Equation 01]), instead, 
would predict -1.47 (1.6a away from ours) and -1.48 (l.ler), in general in reasonable 
agreement with the data. 

From a simple equilibrium equation we can estimate the temperature of the grains, T t 



a- 



as: 



T 9 ~ 



< Quv > 



< Q IR > 16na SB d 2 



40 K: 



(5) 



with < Qnv(jR)-> the bla ckbody-averaged Q-factor in the UV (IR) wavelengths (IDraine 



981 



Draine fc Lee 



((Smith 



19841 ). L s ~ 2 x 10 7 the total stellar luminosity of the Carina nebula 
20061 ) and d ~ 40 pc (1°) the distance of the dust grains from the sources. < Quv > 



is the blackbody averaged absorption efficiency (assuming a blackbody temperature of 
38 x 10 3 K); the typical distance d between the grains and the star has been estimated from 
the angular size of the nebula and its distance. This value, although obtained from a naive 
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model, is in reasonable agreement with our results. 



Equation (JSJ) can be also used to estimate the mean distance of the dust grains from 
the center of the nebula once derived their temperatures. For a dust population of 34.5 K 
we find a distance of about 50 pc; in the presence of a colder dust population of about 18 K 
the distance would be about 200 pc. 



Yonekura et al. 



(120051 ) estimate a mass of the molecular gas of about 3.5 x 10 5 M Q ; 



since the molecular gas in the Carina neb ula is about one t hird of the total mass budget 



(atomic+ molecular gas+ dust) following ISmith fc Brooks! ( 120071 ) we can estimate a total 
mass of about 1O 6 M ; this compared with our estimated mass give s a gas-dust ratio 



of about 100, which i s the canonical value assumed in this nebula ( iMorris et al. 



1999 



Smith fc Gehrz 



2005 



Smith fc Brooks 



2003). 



Our best estimate of the free-free emission, based mainly on the WMAP fluxes, is at 
1.6cr from the Parkes observation. One possibility to explain an excess flux observed in the 
Parkes data would be to suppose the presence of synchrot ron radiation in the nebula, but 



1966 



Jones 



1973). 



this does not fit other long wavelength radio observations (IBeard fc Kerr 
More likely are effects related to the very different angular resolution of the Parkes data, 



2.6 arcmin flHuchtmeier fc Day 



19751 1. with respect to the WMAP observations. With a 
large beam the flux measured by WMAP in the pixels of the box might be overestimated, 
since every peripheral pixel includes the contribution of sources outside the integration box, 
but still inside the fraction of the beam spilling outside the box. The largest WMAP beam, 
0°.81 FWHM, is still smaller than our integration box, so that the flux of the Carina Nebula 
lost in the sidelobes of the beam is negligible. However, if the 0°.8 beam is centered on the 
edge of the integration region, it might pick-up flux from sources outside the integration 
region, where half of the beam solid angle is lying. We have double checked this statement 
with numerical simulations. We have simulated a sky brightness distribution assuming that 
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a single source, 30 arcmin in radius, is located in the center of our integration box. We 
have convolved this brightness with Gaussian beams with the same FWHM as the WMAP 
beams and integrated over our integration box to estimate the measured flux. The result 
is the same within less than 1% (at 94 GHz) (or up to 3% at 23 GHz) as in the case of an 
unsmoothed map, showing that the power lost in the sidelobes is negligible. Then we have 
added a second source, with 30 arcmin diameter and 20% of the maximum brightness of 
Carina, just outside our integration box. In this case we observe an overestimate of the 
measured fluxes which can be, again, from less than 1% (at 94 GHz) to 3% (at 23 GHz) 
of the unsmoothed map. Therefore, the modest resolution of the low frequency channels 
of WMAP does not seem to be the reason of the deficit with respect to the Parkes data. 
Prob ably more important is the fact that the Parkes data do not cover our entire integration 



box ( Smith fc Brooks 



20071 1 . From all this we conclude that a direct comparison is very 
difficult. This issue, however, does not affect the dust properties in the Carina nebula, 
which are the main focus of this paper. 

The Carina nebula SED is well described by a four-parameter fit, although the lack 
of data between 240 and 730 //m cannot exclude the presence of another dust component. 
Our 7-parameter fit gives an hint of this; nevertheless more data in the missing wavelength 
range are necessary to confirm it. 



Our method of averaging over distinct background re gions provides a 



estimation generally in agreement with what was found by 



Smith fc Brooksl (120071 ) 



background 



At wavelengths shorter than 60/im, the emission in the nebula is dominated by 
polycyclic aromatic hydrocarbons (e.g., Rathborne et al. 2002) and warm dust grains 



at temperatures >80K. At wavelengths longer than 60 /xm, the flux 



negligible with respect to the integrated flux of the nebula. 



Cox et al 



r om n Carinae is 



(119951 ) have shown 



that the composite IR SED of r\ Carinae and the Homunculus Nebula is described by two 
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modified Planck spectra i/B(u,Tj) with 5 = 1 and T i= 210, 430 K. At 20 jum, the peak of 
the flux density from the Homunculus nebula is about 10 5 Jy, while about 10 2 Jy comes 



from the 35 K dust population of the more extended Carina Nebula (ISmith fc Brooks 



20071 ): 



at 60 /im the emission of the Homunculus nebula drops to a value of ab out 2 x 10 4 Jy, whi 



that of the 35 K dust increases to about 9 x 10 5 Jy. Among r ecent results 



the ones based on Infrared Space Observatory observations (IMorris et al. 



Smith et al. 



5 



20031 ) 



1999) also confirm 



these findings. The global mass of the Carina nebula is not strongly affected by the mid-IR 
(< 60/xm) portion of the spectrum: the 200 K dust population adds a small contribution to 
the total mass budget which is negligible with respect to the large uncertain coming from 
colder (at 18 or 40 K) dust grains populations. 

Our hint for a second, colder, dust population with a temperature of about 18 K 
would be in good agreement with the analysis of the star-forming regions observed by 
RONAOS, which predic ts a high spectral index (about -2) for very cold regions (11-20) K 



Pupae et al 



2001 



20031 ). The analysis of the millimetric continuous emission performed 



value of the spectral index 



20061 ). However, 



on the SIMBA-SCUBA star- forming regions shows a typica 
equal to -2, assuming dust temperatures higher than 10 K (iHill et al. 
the existence of a cold component cannot be excluded from the available data. It will be 
interesting to improve the SED of the Carina nebula with future observations, both better 
calibrated and at complementary wavelengths, either confirming or disproving hints for the 
existence of the very cold dust component, and its polarization properties. 



2007a 



b) will observe 



The future balloon borne experiment PILOT (IBernard et al. 
the Carina nebula at the two wavelengths of 240 and 550 /im; observations at the first 
wavelength could confirm the flux measured by COBE-DIRBE, observations at the second 
wavelength will fill our missing range. 



The ongoing Planck mission (ITauber et al. 



20101 ) has measured the whole sky at 
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nine wavelengths from 30 to 900 GHz, with unprecedented sensitivity, accurate (< 1%) 
calibration, and with angular r esolution down to 5 arcmin; the SPIRE instrument on 



board of the Herschel satellite (I Griffin et al 



20101 ) will also cover the region at 250, 350, 



and 500 /im. Once processed, these data will certainly solve the issue. The error in the 
background estimates, at wavelengths shorter than 2 mm, are comparable to the difference 
between the SEDs of our one- and two-component best fits and larger than the Planck 
calibration errors; therefore they will not be the limiting factor for discriminating between 
the two models. 



Conclusion 



We have prese nted multi-wavelength o bservations of the Carina nebula. We compared 



our analysis to the ISmith fc Brooks! (120071 ) derivation of the SED. The SED of the nebula 
can be described as the sum of a graybody component, with a power-law emissivity, 
and an optically thin thermal bremsstrahlung. The parameters of the model are well 
constrained by the data: we find e x = (1.10^;1) x 10~ 3 , 7\ = (34.5±?;g) K, a x = -1.37+^1, 
a s = 1113^5 Jy. Thanks to the B98 observations, in the range 730 /zm - 3.3 mm, we have 
been able to derive the spectral index a x rather than assuming it equal to -1. We also 
found, from our multi-component analysis, that the data are consistent with an inverse 
relation between spectral index and temperature. However, the improvement in the \ 2 
is not statistically significant. Our four-parameter results do not exclude the possibility 
that the FIR SED is described by two graybody components, the latter with a colder 
dust temperature and a steeper spectral index than the former. The small number of 
experimental data, compared to the number of free parameters does not allow a robust 
investigation of this hypothesis. This issue will be solved by the Planck/Herschel surveys. 
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Fig. 1. — Wide-field images of the Carina nebula, from far-infrared (top left) to radio (bot- 
tom right) wavelengths . The ranges of the galactic coordinates (in degrees) reproduce the 
Smith fc Brooks! (120071 ) integration box. 
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Fig. 2. — Standard deviation of the mean background contribution, a _F (back ground), to the 
measured flux, F(source), corrected for the background, vs. the wavelength of the different 
data sets; the number (n) labels the different number of regions considered for estimating 
the mean background. The IRAS points are labeled with half-colored symbols. 
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Fig. 3. — Contour plots of the joint probability density and likelihoods for the four parameters 
describing the SED of the Carina nebula. The dotted, dashed and solid regions represent 
confidence intervals for 68.3%, 95.4%, and 99.7% CL, respectively. 



-28 - 




Fig. 4. — FIR to radio SED of the Carina nebula. In the legend we report, from top to 
bottom,: the best fit SED for four- parameter y 2 fit, the sa me for the MCMCMH with four, 



six, and seven parameters, and the 
measured integrated fluxes. 



Smith fc Brooks! (120071 ) SED. The symbols represent the 
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